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This paper is the fourth in a series exploring the physical consequences of the solidity of highly 
viscous liquids. It is argued that the two basic characteristics of a flow event (a jump between two 
energy minima in configuration space) are the local density change and the sum of all particle dis- 
placements. Based on this it is proposed that density fluctuations are described by a time-dependent 
Ginzburg-Landau equation with rates in k-space of the form To + Dk 2 with D 3> Toa 2 where a is 
the average intermolecular distance. The inequality expresses a long-wavelength dominance of the 
dynamics which implies that the Hamiltonian (free energy) may be taken to be ultra local. As 
an illustration of the theory the case with the simplest non-trivial Hamiltonian is solved to second 
order in the Gaussian approximation, where it predicts an asymmetric frequency dependence of the 
isothermal bulk modulus with Debye behavior at low frequencies and an lj _1//2 decay of the loss at 
high frequencies. Finally, a general formalism for the description of viscous liquid dynamics, which 
supplements the density dynamics by including stress fields, a potential energy field, and molecular 
orientational fields, is proposed. 
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I. INTRODUCTION 



Glasses are made by cooling viscous liquids. The liquid relaxation time r increases dramatically upon cooling, and 
the glass transition takes place when r exceeds the inverse cooling rate: r 3> l/\dh\T / dt\. For cooling rates of order 
K/min one speaks about the calorimetric glass transition; glasses, however, may be produced by much higher cooling 
rates like the splat coolings of traditional metallic glasses, or much more slowly as when manufacturing large mirrors 
O \ for astronomical telescopes by cooling gently over months. In all cases, by definition the glassy state has been arrived 
at when the liquid is no more in thermal equilibrium. The glassy state is not unique - it depends on the thermal 
history of the system after it first fell out of equilibrium. In view of this inherent complexity of the glassy state it 
^ ■ appears that a genuine understanding of a glass and its properties can come only from a proper understanding of the 
s ! ' preceding much simpler equilibrium viscous liquid phase. 

Are viscous liquids just quantitatively different from less-viscous liquids or are they, in fact, qualitatively different 
[1-8]? On the one hand, one expects that all liquids obey the Navier-Stokes equation, albeit with a viscosity which 
in some cases is so large that it would take years to pour the liquid out of a cup; this suggests that the difference is 
just quantitative. On the other hand, consider the actual molecular motions in viscous liquids. Going back in time 
at least to Kauzmann's 1948 paper [1], these were always believed to be predominantly vibrational. The physical 
picture is that the molecules are caught in deep potential energy minima, vibrating billions and billions of times before 
rearranging to arrive at another potential energy minimum. Kauzmann described these rare flow events as "jumps 
of molecular units of flow between different positions of equilibrium in the liquid's quasicrystalline lattice" [1]. This 
idea, which was elaborated upon in Goldstein's 1969 paper [4], was confirmed during the last decade in numerous 
computer simulations. In our interpretation the resulting physical picture [9] is: 

Ch , 

q ' Viscous liquid = Solid that flows . (1) 



The implicit statement is: Viscous liquid ^ ordinary liquid. If the vibrations are averaged out, viscous liquid dynamics 
may be identified with the "inherent dynamics" consisting of jumps between the energy minima in the system's 
configuration space [10]. Each minimum has a basin of attraction [4, 5], known as an inherent structure. The concept 
of inherent structures was originally introduced by Stillinger and Weber as a way of thinking about liquids and solids 
in general [11], but the concept seems to be particularly useful for understanding the physics of highly viscous liquids. 

The liquid relaxation time r is related to the (shear) viscosity r\ and the instantaneous shear modulus G M by 
Maxwell's famous relation t = rj/Goo [12]. For "ordinary" less- viscous liquids like ambient water the viscosity is in 
the 10 -3 Pa s range and the instantaneous shear modulus is typically of order 10 9 — 10 10 Pa. Thus the relaxation 
time is barely one picosecond and comparable to a typical molecular vibration time. Viscous liquids approaching the 
calorimetric glass transition, on the other hand, have relaxation times of order 100-1000 seconds. These extremely 
long relaxation times reflect the fact that molecular motion has almost completely seized. The molecules still have 
the thermal velocity distribution, of course, but virtually all motion is vibrational like in a solid. Thus the incoherent 
(single-particle) diffusion constant D s is extremely small, and in the well-known expression for D s in terms of the 
velocity autocorrelation function D s — (v x (0)v x (t)}dt there is a most delicate cancellation of contributions. This 
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fact was emphasized in 1984 by Brawer [7] who pointed out that if one wishes to apply conventional liquid state 
theory to viscous liquids approaching the calorimetric glass transition, the approximations used should be accurate 
to more than 15 digits in order to give reasonable results. 

A further argument for viscous liquids being qualitatively different from the less-viscous liquids dealt with in 
traditional liquid state theory is the following. A liquid is characterized by a number of diffusion constants: The 
heat diffusion constant, the incoherent diffusion constant, the dynamic viscosity v of the Navier-Stokes equation (i.e., 
the transverse momentum diffusion constant: viscosity/density), and the coherent diffusion constant (characterizing 
the long- wavelength decay of density fluctuations). For "ordinary" liquids these diffusion constants are all typically 
within one or two orders of magnitude of 10 _7 m 2 /s. This is easy to understand from kinetic theory, according to which 
the diffusion constant is of order the mean-free path squared over the mean time between collisions of the diffusing 
entity Rough estimates of these quantities are one Angstrom and 0.1 picosecond respectively, resulting in the value 
10~ 7 m 2 /s. For a viscous liquid approaching the calorimetric glass transition, however, the approximate identity of 
diffusion constants breaks down: As the Maxwell relaxation time increases upon cooling from the less- viscous phase, 
the single-particle diffusion constant D s decreases roughly inversely proportionally to viscosity [13]. At the same time 
the dynamic viscosity v increases. Thus the ratio D s /v changes from roughly 1 to a number of order 10~ 30 just above 
the calorimetric glass transition. Such small dimensionlcss numbers are rare in condensed matter physics. Small 
numbers in physics usually signal a qualitative change and a simplification of some kind. It view of this it seems likely 
that the physics of highly viscous liquids are different from - and somehow simpler than - those of "ordinary" liquids. 

Viscous liquids have common features which are independent of the nature of the chemical bonds involved. These 
features [14, 15] may be summarized into three non's: Non-Debye behavior of the main (alpha) relaxation process, 
non-Arrhenius viscosity (or relaxation time, r ~ n because r = ri/G^), and non-linearity of relaxations following 
even relatively small temperature jumps. The last non is probably less basic because in most models non-linearity 
follows from the strong temperature dependence of r, but the non-Debye and non-Arrhenius behaviors constitute 
crucial and defining characteristics of viscous liquids. 

Given the universal features of viscous liquids, an obvious question is: Does the high viscosity - translating into 
(1) - in and of itself make it possible to physically understand and mathematically derive the non's? Since the high 
viscosity is caused by energy barriers large compared to fc^T, this seems to have been Goldstein's view when he in 
1969 wrote: "I am only conjecturing that whatever rigorous theory of kinetics we will someday have, processes limited 
by a high potential barrier will share some common simplifications of approach" [4] . 

The present paper is the fourth in a series (I-III) [16-18] entitled Solidity of viscous liquids attempting to identify 
Goldstein's "common simplifications of approach" by extracting the physical consequences of Eq. (1). The term 
solidity is meant as a concise way of referring to the "solid- like-ness" of viscous liquids. The first paper from our 
group utilizing solidity type arguments preceded the series. This was a joint publication with Olsen and Christensen 
from 1996 [19] which proposed a model for the non-Arrhenius temperature dependence of r (see also Ref. [20]). The 
idea is that the activation energy for a flow event is the work done in shoving aside the surroundings in order to - in a 
brief moment - create the extra space which by assumption is needed for the molecules to rearrange. According to this 
"shoving" model the activation energy is mainly shear elastic energy located in the surroundings of the rearranging 
molecules. Consequently, the activation energy is proportional to the instantaneous shear modulus, a quantity which 
is usually much more temperature dependent in viscous liquids than in crystals and glasses (or in "ordinary" liquids) . 
The shoving model appears to work well for molecular liquids [19], but it is too early to tell whether the model is 
generally applicable [15]. 

Paper (I) introduced the concept of a solidity length I that determines the length scale below which a viscous liquid 
for all purposes behaves as a solid, albeit one that flows. In terms of the average intermolecular distance a (defined 
by writing the volume per molecules as a 3 ), the alpha relaxation time t, and the high-frequency sound velocity c, the 
solidity length is given by 

I 4 = a 3 rc. (2) 

This expression was derived by noting that a flow event is followed by the emission of a spherical sound wave; I is 
determined by requiring that elastic equilibrium is just about established throughout a sphere with radius I before 
the next flow event inside the sphere typically takes place [16]. At the calorimetric glass transition the solidity length 
is several thousand Angstrom. 

In (II) we discussed anisotropic flow events and defined a parameter characterizing the anisotropy; it now appears, 
however, that this parameter is of only minor significance - see below. The origin of the non-Debye alpha linear 
response functions was the subject of paper (III) where it was argued that a long-time-tail mechanism may explain 
what appears to be a generic property of the alpha process as observed, e.g., by dielectric relaxation [21]: At high 
frequencies the loss follows an cj -1 / 2 decay once the effects of beta processes are minimized by going to sufficiently low 
temperatures (but still in the equilibrium liquid phase). It was further argued that the coherent diffusion constant 
D is much larger than the incoherent diffusion constant D s , implying that D 3> a 2 jr. The latter inequality, which is 
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essential for the long-time-tail mechanism to work for a range of times shorter than r, i.e., above the alpha loss peak 
frequency, reflects a long- wavelength dominance of the dynamics. 

Long-time tails usually derive from a diffusion equation which implies decay rates oc k 2 for the k-wavevector 
component of the relevant conserved variable [22, 23]. An explicit realization of the long-time-tail scheme of paper 
(III) was given subseqently in a simple model for dielectric relaxation [24] . This model is based on a conserved scalar 
field, the density, and a non-conserved vector field, the dipole density, with the simplest possible interaction term. 
When the model is solved in the Gaussian approximation to second order in the interaction strength between the two 
fields, one finds for the frequency-dependent dielectric constant 

1 1 c 
e(u) oc - + -= + : — . (3 

VI + «WT V2 + VI +ILOT 1+UVT 

This expression reproduces the generic [21] asymptotic behaviors of the alpha process: — e" awat low frequencies 
and — e" cx 1/\/uj at high frequencies. 

The present paper elaborates on the point made briefly in (III) and Ref. [24] that the density "dispersion relation" 
(the decay rate) realistically should include an additive constant and be of the form T + Dk 2 where T = 1/r. We 
arrive at this result by first identifying the most important characteristics of a flow event. A new argument is given for 
the long- wavelength-dominance inequality D 3> Toa 2 . Then a model for density fluctuations is proposed and solved 
in the simplest approximation, predicting that the loss peak of the frequency-dependent isothermal bulk modulus is 
asymmetric towards the high-frequency side where the loss decays as oc a;" 1 / 2 . Finally, we propose principles for the 
general description of viscous liquid dynamics. 



II. BASIC CHARACTERISTICS OF A FLOW EVENT 



A. A first approach 

The arguments of this paper refer to length scales below the solidity length I. For liquids approaching the calorimctric 
glass transition I is so large that there are almost four orders of magnitude of length scales between I and the 
intermolecular distance a. We argued above that in view of Eq. (1) viscous liquid dynamics are basically to be 
identified to the "inherent dynamics" [10] consisting of jumps between potential energy minima (inherent structures). 
This is what Goldstein envisaged in his emphasis on the importance of potential energy barriers much larger fc#T 
[4]. The fact that these barriers are large is directly responsible for the extremely large low- temperature viscosity, 
of course, but it also implies that there is a clear separation between the molecular vibrations and the much slower 
changes of the configurational degrees of freedom associated with jumps between the potential energy minima. In the 
resulting picture the vibrations are regarded as uncorrelated to the inherent dynamics, and only the latter contribute 
to autocorrelation-function variations on the time scale of the alpha relaxation. 

A flow event is a jump from one potential energy minimum to another. Leaving aside the interesting question why 
the alpha-process activation energy increases upon cooling, and whether this as in the shoving model [19] reflects the 
fact that the high-frequency elastic constants increase upon cooling, we shall nevertheless use shoving-model type 
arguments below. This is done by comparing the situation before and after a flow event, and not as in the shoving 
model before and at the barrier maximum of the flow event. 

It seems compelling that flow events are localized in space; this is also what is observed in the numerous computer 
simulations that now have been published. Justified by Eq. (1) we utilize arguments from the theory of solid elasticity 
[25] by regarding the liquid as an isotropic solid in which the flow event takes place. As a simple model, suppose 
that the flow event is radially symmetric. If it takes place at r = 0, the radial displacement field induced in the 
surroundings may be written u(r) = u r (r)v/r where r = |r|. To determine u r (r) we combine the equation of elastic 
equilibrium after the flow event - zero divergence of the stress tensor: <9i(7ifc = - with the stress-strain relation 
Uik = Kuu6ik + 2G(uik — uuSik/3), where K and G are the bulk and shear moduli and Uik = (diitfc + dkUi)/2 is the 
strain tensor. This leads to 

(A" + -)V(V ■ u) + GV 2 u = 0. (4) 

Because the displacement is radial one has V x u = which, via the vector identity V x (V x u) = V(V • u) — V 2 u, 
implies that V(V • u) = V 2 u. When substituted into Eq. (4) this leads to V(V • u) = 0, or V • u = C\. The radial 
displacement is found by solving the equationV • u = r~ 2 d r {r 2 u r ) — C\, leading to u r — C^r -2 + C\r/3. Since 
the latter term diverges as r — > oo, we must have C\ = 0. In conclusion, the displacement field is analogous to the 
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Coulomb electric field of a point charge: 



r 

u oc — . 



(5) 



The relative density change is — V • u [25] which is zero except at r = 0. Thus in this macroscopic and radially 
symmetric description there is no density change in the surroundings of the flow event. A decrease of density at 
the flow event center induces a positive radial expansion which results in the same particle flux through all spheres 
centred at the flow event. This means that for any spherical shell surrounding the flow event, thick or thin, just as 
many molecules enter the shell from the inside as leave it on the outside. In effect, a flow event corresponds to a 
three-dimensional version of Hilbert's hotel, the infinite hotel which even when totally occupied - makes room for 
an extra guest by asking all guests to move to one higher room number. 

The above analysis is oversimplified because of the assumption of spherical symmetry. Nevertheless, the result that 
the displacement field far from the flow event center varies as 1/r 2 is correct and general. To see this, note that the 
effect of one flow event on its surroundings may be mimicked by first imagining a tiny sphere surrounding the flow 
event center. If the molecules inside the sphere are removed, the effect of the change of positions of the molecules 
before and after the flow event may be reproduced by external forces acting suitably on the surface of the sphere - 
forces that must sum to zero. An external force acting on a point in an elastic medium introduces a momentum flux 
spreading to infinity; since the stress tensor a is the momentum flux density, we conclude that cr oc 1/r 2 as r ^ oo. 
The stress tensor is given by first order derivatives of the displacement vector field; thus one expects that |u| cx 1/r. 
For forces summing to zero, however, we get |u| cx 1/r 2 , just as the potential from an electric dipole varies as 1/r 2 
whereas it varies as 1/r from a point charge. As shown elsewhere [26], even in the most anisotropic case more than 
90% of the elastic energy in the far-field surroundings is shear elastic energy, so the assumption that there arc no 
density changes in the surroundings is a good approximation for describing the long-ranged effects of a flow event. 

Flow events may be regarded as instantaneous on length scales below the solidity length. Numbering the flow 
events consecutively after the time they take place, i M , if r M is the center of the p'th flow event and the number 
measures its magnitude, the above considerations translate into the following dynamic equation for the density p(r, t) 
in a coarse-grained description: 



Equation (6) is not inconsistent with particle conservation, of course. Nevertheless, a field obeying Eq. (6) does 
have the appearance of not being conserved, because density changes at one point in space do not affect the density 
elsewhere. A field of non-interacting spins fluctuating randomly in time would be described by a similar time evolution 
equation, and there is nothing conserved by such a spin field. To summarize, the "Hilbert's hotel effect" deriving 
from solidity is not inconsistent with the fact that the particle number is obviously and trivially conserved, but it 
implies that density acts as if it is a non-conserved field. Although we show below that Eq. (6) is too simple to 
reflect all relevant features of density fluctuations, this conclusion remains valid. The interpretation of the flow event 
magnitude & M is found by integrating Eq. (6) over a brief period of time including only the flow event taking place 
at i M . This leads to a density change equal to b^S(r — r M ), so the number of particles inside any volume which 
includes the flow event center r M changes by precisely & M . Note that via the continuity equation p + V • J = and the 
identity V • (r/r 3 ) = 47r<5(r), Eq. (6) corresponds to the following expression for the particle current density in the 
coarse-grained description 



To study flow events in more detail we consider the induced density changes by going to k-spacc. If the liquid 
consists of N molecules in volume V each with position r j , the variable pk defined by 




(6) 




(7) 



B. 



A more detailed treatment 




(8) 



Normalizing in this way is convenient because it makes the fluctuations independent of sample size in the V — > oo 
limit where the static structure factor is given by S(k) — (|pk| 2 ) [27]. 
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As usual, the k- vectors are restricted to values compatible with periodic boundary conditions, i.e., having an integer 
number of periods in the volume V in all three axis directions. We are particularly interested in small k-vectors. If 
the molecular displacements induced by a single flow event are denoted by Ar,- , by a first order Taylor expansion the 
change of pk for small k is given by 

^ = ^E elkr ^ k - A ^- ( 9 ) 



Since particle displacement are unlikely to be much larger than the intermolecular distance a, Eq. (9) applies whenever 
ka <C 1. As discussed in (III) and Ref. [28], for realistic viscous liquid samples momentum is not conserved because 
the transverse momentum diffusion constant (the kinematic viscosity v — f]/(p)) is so large that on the time scale of 
the alpha relaxation momentum is unavoidably exchanged between the sample holder and the liquid. Thus momentum 
conservation is irrelevant for viscous liquid dynamics just as it is, e.g., for the description of point defect diffusion 
in crystals, and the sum of all particle displacements induced by a flow event AR = 5^ . Arj is generally non-zero. 
In our present view the vector AR is a more important measure of flow event anisotropy than the "quadrupolar" 
parameter introduced in (II), where it was implicitly assumed that AR = 0. 

If the /z'th flow event is centered at r M , because the largest displacements take place close to r M , for small k it is 
tempting to argue (III) that it is a good approximation to replace all exponentials by exp(zk • r M ), leading to 

ik-r 

S Plc = < —=Lik ■ AR. (10) 
V N 

This implies that |(5/?k| 2 oc k 2 for small k. Equation (6), on the other hand, which we argued provides a good coarse- 
grained description, implies that |(5/?k| 2 is constant for small k. Which is right? In fact, both are partially correct as 
we now proceed to show. 

The starting point is Eq. (9) (still assuming small k). The calculation leading to Eq. (10) is not quite correct, 
however. This is because, although displacements far from the flow event are small (oc 1/r 2 ), there are many molecules 
far away (oc r 2 ) and their contributions cannot be ignored. A more detailed analysis proceeds as follows. Suppose a 
flow event of magnitude 6 is located at r = 0. As a reasonable first approximation we use Eq. (7) which implies that 
the total particle flux due to this flow event through the perpendicular area dA at the distance r from is equal to 
(— &/47r)(cL4/r 2 ). Identifying this flux with {p)udA, where u is the particle displacement, leads to u = — (b / Anr 2 {p}) . 
Thus the displacement of the j'th molecule, Ar^, is given by Ar, = — (b/4n(p))(rj/r 3 ). In terms of the density 
p(r) = 5(r — Vj) Eq. (9) thus becomes 

. -b f p(r) ik ■ r ik 

5p k = = / ^^-e 4kr dr . 11 

4nVN Jv (P) r 3 

If the flow event magnitude b is uncorrelated to other quantities, the ensemble average of the absolute square is given 

At first sight this expression appears to confirm that (|<5pk| 2 ) oc k 2 , but that is not correct: As shown in the Appendix, 
because the term (p(r)p(r')) for |r — r'| — > oo becomes constant, one has 



k' 



For small k this implies that 

where C ~ (l/N) £ k , S(k')/[(k'a) 2 } ~ 1. 



(\Sp k \ 2 ) oc 1 + C(ka) 2 , (14) 



III. DENSITY FLUCTUATIONS 
A. General framework 



Based on the above considerations we now seek a model for density fluctuations in equilibrium viscous liquids. 
To simplify matters it is assumed that density is the only relevant variable, although eventually other variables 
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should definitely be included in the description (Sec. IV). The basic assumption is that the density obeys a standard 
Langcvin equation [29]. For discrete degrees of freedom Qi, ...,<3„ the Langevin equation starts from a Hamiltonian 
H(Qi, ...Q n ), and the dynamics are given by the equations Qi = —Tid((3H)/dQi + £j(t) where [3 = l/fegT and £j(i) 
is a Gaussian white noise term obeying (£,i(t)tij(t')) = 2TiSijS(t — t'). These equations give the correct canonical 
equilibrium probability for the average occupation in configuration space [29]. 

In the present case the Hamiltonian is the free energy written as a functional of the density field. In terms of the 
complex density field variables pk the Langevin equation looks as follows with a complex noise term: 

Ak = -r fc ^^ + &(i). (15) 

The noise term obeys = £-k(i) and (£k(i)£k' (*')) = ^kSk,k'5(t - t'); because p k = p_k Eq. (15) is equivalent 

to two independent real Langevin equations, one for the real part of p k and one for its imaginary part. To determine 
the k-dependence of the rate we note that when Eq. (15) is integrated over a short time interval, the noise term 
dominates. Thus if the short time interval is At, the magnitude of the change Apk is given by (|Apk| 2 ) = 2TkAt. On 
the other hand, because flow events are uncorrelated over short time spans, from Eq. (14) one finds that (|Apk| 2 ) oc 
At(l + C(fca) 2 ) for small k. In conclusion, for small k the rate 1^ is of the form [31] 

r fe = r + Dk 2 . (16) 

We shall assume that this expression applies for all k. A conserved field is characterized by Tk oc k 2 for k — > [30], 
so Eq. (16) expresses the fact that density has the appearence of a non-conserved field. According to the calculation 
of last section C ~ 1, which implies that D ~ F a 2 . This calculation, however, assumes much more symmetry than 
realistically may be expected. Thus Eq. (5) is based on the macroscopic elasticity theory describing an isotropic and 
homogeneous solid. Surely a viscous liquid is solid like, but it is neither homogeneous nor isotropic on the molecular 
scale, and although the fact that the displacements far from a flow event vary as 1/r 2 remains valid, one wouldn't 
expect Eq. (5) to be accurate. Violations of isotropy and homogeneity easily lead to a much larger flow event induced 
AR than found in the calculations leading to Eq. (14) (similarly, local correlations of flow events as seen, e.g., in 
the observations of strings of flow events in some computer simulations [10, 32], also severely violate isotropy). In 
conclusion, since the C term of Eq. (14) is determined by the magnitude of AR which is a measure of the flow event 
anisotropy, it appears likely that C > 1, or cquivalently 

£>>r a 2 . (17) 

We shall henceforth assume that this inequality is obeyed. In (III) Eq. (17) was justified by arguing that D is much 
larger than the incoherent diffusion constant. In neither case, however, has Eq. (17) been rigorously proved and it 
remains an assumption which is justified by physical arguments. 

In principle, all k-vectors consistent with periodic boundary conditions are allowed, but since there are just a 
finite number of molecules, one should only include the N smallest k-vectors - for larger k-vectors the pk's become 
redundant. This means that there is an implicit cut-off in k-space, k c , which is easily shown to be given by k c a ~ 1. In 
view of this, Eq. (17) implies that there is a range of allowed k-vectors where the Dk 2 term of Eq. (16) dominates over 
the Tq term. As shown below and in Ref. [24] this makes it possible to understand the generic uj^ 1 ^ 2 high-frequency 
decay of the alpha loss as a consequence of a long-time tail mechanism operating over a range of times short compared 
to the alpha relaxation time. 

The inequality Eq. (17) expresses a long-wavelength dominance of the dynamics. This assumption, which a posteriori 
justifies our focus on small k vectors, makes it possible to simplify the Hamiltonian considerably. A number of 
scattering experiments have looked for diverging length scales as the glass transition is approached, but found none. 
The consensus is that viscous liquids have no long-ranged static (i.e., equal time) density correlations. Thus, if the 
dynamics are dominated by the long-wavelength behavior, a model with no equal-time spatial density correlations 
should suffice. In field-theory terms this means that the Hamiltonian may be assumed to be "ultra local," i.e., 
without the usual gradient term (Vp) 2 or other terms coupling fields at different points in space. Although this 
is rather unusual from a general field-theory perspective, note that a simple example of an ultra-local field theory 
is the free energy functional for an ideal gas. More generally, the Ramakrishnan-Yussouff density functional [33] 
becomes ultra local in descriptions which are coarse-grained on length scales beyond the correlation length of the 
direct correlation function. 

Thus we assume that the dimensionless Hamiltonian functional, when scaled by the inverse temperature (3, is of 
the form 



0H = (p) / dr 



v 



1 (p{v)-(p)\ 2 , A fp(r)-( P y 3 



2A V ( P ) J 3 V (p) 



(18) 
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The third order term is crucial for the results derived below; to ensure stability, however, there must be further 
higher order even terms. To transform this expression into k-space, we first note that Eq. (8) implies that p(r) = 
(V~N /V)J2kPkexp(— ik • r), which in turn implies that (p(r) - (p) )/{p) = J^k^o P k cx p(~ *k • r)/\/~N. When this is 
substituted into Eq. (18), the dimensionless Hamiltonian becomes 

& H = jA.J2 pkp - k + X^ k ^ k '^- k - k ' + ( 19 ) 

where it is implicitly understood that no terms with k = appear (because a system with fixed volume is considered, 
for k = /?k is not a dynamic degree of freedom). 

Once the dynamics have been specified one can calculate the density autocorrelation function as a function of time, 
the Laplace transform of which gives the dynamic structure factor. To the best of the author's knowledge, there are 
no data for this quantity for highly viscous liquids. But by the fluctuation-dissipation theorem the k — > limit of 
the density autocorrelation function determines the macroscopic frequency-dependent isothermal bulk modulus. It is 
convenient to introduce the notation 

*k(i) = (pk(O)p-k(i)) • (20) 

Obviously, &k(t) = 3>_k(f) by time-reversal and parity invariance. It is possible to establish an exact equation for 
®k(t) by use of the following general theorem [18, 24]: If Qi obeys a Langevin equation of the form (no sums over 
i) Q, = -Tidi{l3H) + &(*), one has d 2 /dt 2 (Q t (0)Q l (t)) = T 2 (d^i?) (())<% (/3i?)(t)). Substituting Eq. (19) into this 
identity leads to 

. r^(^ + ^S A+k ,(0),_ l ,(0 ) + ...)(4<fi + -^S^--('K»W + ^)) ■ (2D 

Approximations are needed in order to proceed. The simplest approximation is the Gaussian approximation which 
leads to non-linear self-consistent equations. 



B. The Gaussian approximation 

The simplest nontrivial case is when one ignores the higher-order terms of Eq. (19), a procedure which is jus- 
tified when these terms are so small that they do not significantly influence the density autocorrrelation func- 
tion. Using the well-known fact that for variables with zero mean distributed according to a Gaussian one has 
(X1X2X3X4) = {xiX2){x 3 ,x i ) + (xiX3)(x2X4) + {xix i ){x2X^} whereas averages of odd order are zero, Eq. (21) in the 
Gaussian approximation becomes 

^*k(t) = II + 2^X> k+k ,(i)d>_ k ,(t) + • (22) 

We keep only the terms relevant for the below calculation; as mentioned the inclusion of, e.g., a small fourth order 
term in the Hamiltonian may be ignored for calculating the density autocorrelation function - its inclusion would 
simply lead to a minor renomalization of A and a small <J> 3 -type term. 

The k — > limit of $k(i) may be determined analytically to second order in A by proceding as follows. To zeroth 
order in A the Hamiltionan is quadratic, implying for the equal time average that (pkP-k) = A. Thus to zeroth order 
one has $ k °^(i) = Aexp(— T^t/A), which is substituted into the perturbing term of Eq. (22). Because of Eq. (17) the 
k-sum may be evaluated by extending the k-intcgration to infinity (recall that (p) = a~ 3 ): 



^k e -2(T a +Dk' 2 )t/A 



f 

2X1 ^ - e - 2r °</^ r dk > 47 rk' 2 e - 2Dk ' 2t / A 
Jo 



{2nf(p) 

2X 2 A 2 , (2_pty^Vl (23 ) 

\ 2 A 7 /2 

1 _ t -3/2 e -2T t/A _ 



8y/2(nD) 3 / 2 (p) 
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If the k — > limit of $u(t) is denoted by $(t), we thus arrive at the following equation for <&(t) in terms of the 
dimensionless time t = T t/A 

*(t) = + A?- 3 / 2 e- 2 *, (24) 
where A = A 2 A 4 (r /7rD) 3 / 2 /(8v / 2(/o))- The general solution of this differential equation obeying $(t — » oo) = is 

sinh(f - t)e~ 2i 't'- 3 / 2 dt' + C e"* . (25) 

This result applies to second order in A and assumes that the first term is small, but at this stage it is not possible 
to realistically estimate the relative weights of the two terms. 

According to the fluctuation-dissipation theorem, if v is a large subvolume of V, Av(t) = v (t) — v (0), and L(s) is 
the Laplace transform of (Av 2 (t)} evaluated at the complex variable s = iu>, the frequency-dependent isothermal bulk 
modulus Kt(uj) is given by 

tfrMoc-^. (26) 

Besides the relaxational density fluctuations described by Eq. (15) there are always the fast vibrational density 
fluctuations. In any reasonable model these two types of fluctuations are uncorrelated. Consequently, if Av r (t) is the 
relaxational volume change over time t and Av v (t) the vibrational analogue, one has (Av 2 (t)} = (Av 2 (t)) + (Av 2 (t)). 
On time scales much longer than phonon times the vibrational volume mean-square displacement is independent of 
time and one may write (Av 2 (t)) = C + (Av 2 (t)). Thus L(s) becomes C/s plus the Laplace transform of (Av 2 (t)), 
a quantity that is proportional to $(0) — &(t)- Equation (25) implies that ^(O) — $(t) cx \fi at short times (i.e., 
whenever t <C 1 or equivalcntly t <C r where t = A/To), whereas for J > t $(£) goes exponentially to zero. These 
asymptotic behaviors imply that Kt(ui) oc 1 — C\/^/iurr for ljt 3> 1 and Kt{u) cx 1 + Ciiiwr) for lot <C 1. For 
the imaginary part of the frequency-dependent bulk modulus (the loss) one finds K!^(oS) cx lj^ 1 / 2 for lot 3> 1 and 
K'JpiiS) oc cj for lot <C 1. The model thus predicts bulk modulus loss peaks which are asymmetric towards the high- 
frequency side. This is what is always observed for the dielectric and shear mechanical loss peaks, but unfortunately 
there are only few measurements of the frequency-dependent adiabatic bulk modulus (and none of the isothermal 
frequency-dependent bulk modulus). The only published adiabatic measurement known to this author (on glycerol 
[34]) is not inconsistent with these predictions - it is described by a stretched exponential relaxation function with 
exponent 0.43 which implies an asymmetric loss peak fairly close to that predicted here. 



IV. GENERAL DESCRIPTION OF VISCOUS LIQUID DYNAMICS: A PROPOSAL 

To simplify matters as much as possible this paper focused on density as the sole relevant configurational variable. 
In a recent work [24] we discussed the case where the relevant fields are the density and the dipole density fields. This 
was suggested as a simple model for dielectric relaxation. In that model the density and dipole density fields couple 
to each other by a third order term in the Hamiltonian (model C of Rcf . [30] ) . In most cases the molecules do not 
have continuous rotational symmetry around one axis, however, and their orientations should properly be described 
by more general variables CI € SO (3) representing the Eulerian angles. The stress tensor is another relevant field 
which should be included in the description. Finally, the potential energy density should also be included. 

Inspired by the model for density fluctuations discussed in this paper and the above-mentioned model for dielectric 
relaxation [24], we propose the following general recipe for modelling viscous liquid dynamics: 

1. The relevant degrees of freedom are fields 4>^\r), (r) defined as: a) the densities of the different types 
of molecules, b) the densities of the different molecules' orientational variables reflecting their symmetry, c) the 
six stress tensor fields, and d) the potential energy density. 

2. The Hamiltonian H (the free energy) is ultra local; H consists of invariant (i.e., scalar) terms up to some even 
order. 

3. For each field the dynamics are described by a Langevin equation, 



where $ (t) is a standard Gaussian white noise term. 
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4. For each density field the Langevin equation coefficients are given by = T { 3) + D^k 2 where > if } a 2 , 
for all remaining fields the rates are k- independent: — TqK 

V. SUMMARY AND DISCUSSION 

At low temperatures viscous liquid dynamics may be identified with the inherent dynamics taking the system from 
one potential energy minimum to another. Each such jump is referred to as a flow event. The basic characteristics of 
a flow event are: 1) Its position; 2) Its time; 3) The induced density change at the position of the flow event quantified 
by the flow event "magnitude" b of Eq. (6); 4) The total displacement of all molecule positions induced by the flow 
event, AR. 

The density dynamics discussed in this paper are described by a time-dependent Ginzburg-Landau equation [30] 
with rates in k space of the form r + Dk 2 . The first term, which is a consequence of the solidity of viscous liquids, 
reflects the fact that the density may change locally without changing in the surroundings, a result which implies that 
density has the appearance of a non-conserved field variable. The second term, Dk 2 , is the standard diffusion term. It 
was argued that the disorder of the solid-like liquid and its lack of isotropy most likely result in anisotropic flow events. 
This is mathematically reflected in the inequality D 3> Toa 2 which implies an important long-wavelength dominance 
of the dynamics. Since experimentally there are no long-ranged static density correlations, the long-wavelength 
dominance of the dynamics makes it realistic to assume that the Hamiltonian is ultra local. 

Glarum in 1960 suggested that relaxation takes place via defect motion [35, 36]. In his words "molecules do not 
relax independently of one another, and the motion of a particular molecule depends to some degree on that of its 
neighbors... This is because the reorientation of a molecule is far more likely immediately after one of its neighbors 
has relaxed than it is at an arbitrary time." The defect motion was described by a diffusion equation. Glarum's 
model gave a novel mechanism for explaining dielectric loss peaks which are asymmetric towards the high-frequency 
side. Inspired by this work Anderson and Ullman in 1967 generally considered the effect of a fluctuating environment 
on molecular relaxation rates [37]. These authors showed that fluctuations fast compared to the alpha relaxation rate 
lead to a symmetric almost Debye response and slow fluctuations to a symmetric, but broad loss peak - two results 
which were known already from Kauzmann's 1942 paper on dielectric relaxation [38] - whereas if the environment 
fluctuates on the alpha time scale, the loss peak becomes asymmetric towards the high-frequency side. Montrose and 
Litovitz in 1970 [39] discussed a model involving an order parameter with both a diffusive and a decaying term in its 
dynamics. The decay term derives from the fact that "the structure can change by a simple rate process," whereas 
if for instance the order parameter is a function of the number of holes, the diffusive term corresponds to the "many 
small steps arising from the rapid jiggling of molecules." In their review from 1972 of the mechanical response of 
viscous liquid [40] Dexter and Matheson summarizes these three papers by the remark: "Thus, the physical basis of 
the theories of Glarum, Anderson and Ullman, and Montrose and Litovitz is similar: the molecular environment is 
assumed to change as a result of spontaneous molecular motion and small diffusional motions." The similarity to 
the present work is clear; here, however, the "small diffusional motions" occur in the surroundings of a flow event, 
and may be thought of as consequences of the "spontaneous molecular motion." After these early works, also in the 
spirit of the present paper Zwanzig [41] and MacPhail and Kivelson [42] in the 1980's explored the possibility that 
long-time tail mechanisms are relevant for understanding viscous liquiid dynamics. We mention these works in order 
to emphasize that the ideas of (III) and this paper have close analogs in several papers published a long time ago. 

In this paper we regarded flow events as taking the system instantaneously from one inherent structure to another. 
Although this approximation applies only on length scales below the solidity length, which is of order 5,000 A around 
the calorimctric glass transition, it is possible that the proposed density dynamics apply also at macroscopic length 
scales. The problem of properly linking the behavior on length scales below the solidity length to those above needs 
further consideration, though. 
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APPENDIX 



By writing (p(r)p(r')) = (p) 2 + (Ap(r)Ap(r')) where (p) = N/V and Ap = p — (p), Eq. (12) is split into two terms: 



(I^Pkl 2 ) = 



(b 2 ) 
16tt 2 A 



h{k) + h{k) 



where 



and 



h 



(*) = / / 
Jv Jv 



k • r k • r' 

v jv r3 r ' 3 



k • r 



lkr dr 



(A.l) 



(A.2) 



(A.3) 



The latter integral may be expressed in terms of the static structure factor, which for k ^ is defined by 
S(k) = (pk/O-k) = (|pk| 2 ): First note that the definition of /?k (Eq. (8)) implies that for k ^ one has 
v^pk = J v Ap(r) exp(ik • r)dr, which implies that S'(fe) = 1/iV / y J v (Ap(r)Ap(r')} exp(ik • (r - r'))drdr' = 
V/iV J y (Ap(O)Ap(r)) exp(ik • r)tfr. Inversion of this Fourier integral allows one to write the density autocorrela- 
tion function as a sum over k vectors consistent with periodic boundary conditions: 



(Ap(r)Ap(r')) _ 1_ v V k'.( r - r ') 
When this is substituted into Eq. (A.3) we get 

Thus in terms of the integral 



(A.4) 



j(k,n) ee / 



we have 



(l^kl 2 ) = 



16tt 2 A 



|/(k,k)| 2 + 1^5(fc')|/(k,k + k')| 



(A.5) 



(A.6) 



(A.7) 



To evaluate I(k, SI) we note that, if the z-axis is along the SI vector, k is in the xz-plane, and 9o is the angle between 
SI and k, we have SI — (0,0,0) and k = (k sin# , 0, k cos 6> ). In spherical coordinates we thus get 



/■7T />00 f- 

J(k,n) = / d9sin9 / drr 2 

Jo Jo Jo 



This reduces to 



27r /cr(sin# cos</>sin# + cos#cos# ) jOr-cosS 



iOr cos # 



I(k, SI) = 2nkcos9Q I d9sin9cos9 / dre 



The radial integral is evaluated by assuming an implicit convergence term, lim^oo exp(— r/a), leading to 

1 

I cos 9) 



r i 

ICk, SI) = 2nkcos9o / d9 sin 9 cos 9- — — — = Atti — cos#o, 

J (-iSlcr-'-- 



k 



or 



k • SI 

j(k,n) = ^-^r 



(A.i 



(A.9) 



(A.10) 



(All) 
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In conclusion, Eq. (A. 7) becomes 



{\Spk\ ) = — 



l SS (,y k - (k+k ' ; 

k' 



(k + k') 2 



(A.12) 
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